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Segments with the amino acid sequence EKAYLRT appear in natural occurring proteins both in 
a-helices and /3-sheets. For this reason, we have use this peptide to study how secondary structure 
formation in proteins depends on the local environment. Our data rely on multicanonical Monte 
Carlo simulations where the interactions among all atoms are taken into account. Results in gas 
phase are compared with that in an implicit solvent. We find that both in gas phase and solvated 
EKAYLRT forms an a-helix when not interacting with other molecules. However, in the vicinity of 
a /3-strand, the peptide forms a /3-strand. Because of this change in secondary structure our peptide 
may provide a simple model for the a — > /9 transition that is supposedly related to the outbreak of 
Prion diseases and similar illnesses. 

I. INTRODUCTION 

Despite considerable progress over the last decade the problem of predicting the biological active structure of a 
protein solely from the sequence of amino acids has remained a formidable problem. More successful have been 
attempts to predict only the secondary structure. Given the protein sequence it is today possible to determine the 
distribution and location of a-helices and /3-sheets with up to 90% probability. This high success rate indicates a close 
relation between sequence information and secondary structure. However, two observations indicate that this relation 
is not a simple one. First, certain sequences can form either a-helices or /3-sheets [1]. The most prominent example is 
the 11-residue Chameleon peptide [2] that folds as an a— helix when replacing residue 22-32 of the primary sequence 
of the IgG-binding domain of protein G (57 amino acids), but as a /3-strand when inserted instead of residues 42-52. 
Secondly, it has become clear over the last years that miss- folding of proteins, often involving formation of /3-sheets 
instead of a-helices, and subsequent aggregation is the cause of various illnesses including Alzheimer's disease, BSE 
and other Prion diseases. Hence, it is important to understand in detail how secondary structure formation and its 
role in the folding process depends on the intrinsic properties of the protein and the interaction with the surrounding 
environment. 

In order to study these questions we have simulated a peptide whose sequence of amino acids EKAYLRT (glutamine 
- lysine - alanine - tyrosine - leucine - arginine - threonine) appears in natural occurring proteins with significant 
frequency at positions of both a-helices and /3-sheets. Our present work differs therefore from previous investigations 
where we have focused on helix- formation and folding in homopolymers and artificial peptides [3-8]. Unlike these 
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molecules that have a strong intrinsic tendency to form one specific kind of secondary structure elements (a- helices), 
EKAYLRT allows one to research the selection of either helix or sheet, or the transition between these two secondary 
structures, as a function of external factors. 

Our work differs from similar approaches [9,10] in that we study not minimal models but simulate detailed repre- 
sentations of our peptides where the interactions between all atoms are taken into account. EKAYLRT is simulated 
both in gas phase and with an implicit solvent. Quantities such as energy, specific heat, sheetness and helicity are 
calculated as functions of temperature. We find that both the solvated molecule and EKAYLRT in gas phase form 
an a-helix when not interacting with other molecules. However, in the vicinity of a /3-sheet the peptide prefers also 
to form strand. Because of the resulting "auto-catalytic" property our peptide may therefore provide a simple model 
for the a ^ (3 transition and the resulting aggregation process in some proteins that supposedly is related to the 
outbreak to neurological diseases such as Alzheimer's and the Prion diseases. 



Our aim here is to research how secondary structure formation and its role in the folding process depend on either 
the intrinsic properties of a protein or its interaction with the surrounding environment. For this purpose, wc have 
considered detailed, all-atom representations of peptides that are based on the sequence of amino acids EKAYLRT. 
To be more specific, the peptide NH2-EKAYLRT-COOH is studied both as an isolated molecule and interacting with 
another EKAYLRT peptide that is held in a /3-strand conformation. Since our program package SMMP [12] in its 
current version allows only the simulation of single peptides we have modeled the latter case by considering the peptide 
NH2-EKAYLRT-GGGG-EKAYLRT-COOH, with the C-terminal EKAYLRT residues kept as a /3-strand. The four 
glycine residues form a flexible chain that hold the two peptides together but allows their relative positions to vary. 
The underlying assumption is that the interaction between the two EKAYLRT chains is the dominant term and their 
interaction with the glycine residues can be neglected. 

The intra-molecular interactions are described by a standard force field, ECEPP/3, [11] (as implemented in the 
program package SMMP [12]) and are given by: 



II. METHODS 



Eecepp/3 = + Eydw + Ehb + Et, 



(1) 




(2) 
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KdW = -TT - , (3) 



I „12 „6 I 

Ehb = 2^ -12-71^ 



(4) 

Etor = 5Z C/; (1 ± cos(naO) • (5) 

I 

Here, (in A) is the distance between the atoms i and j, and xi is the Z-th torsion angle. The peptide bond angles 
are set to their common value ui = 180°. We further assume for the electrostatic permittivity in the protein interior 
£ = 2 (its common value in ECEPP simulations). 

Simulations of our peptide EKAYLRT in gas phase are compared with such where the interaction of the peptide 
with surrounding water is approximated by an implicit solvent [13]: 

E = EecEPP/3 + Eaolv with Esolv = ^ f i^i • (6) 

i 

Here, Egdy is the solvation energy and thought to be proportional to the solvent accessible surface area Ai of the ith 
atom. The parameters cr, are experimentally determined proportionality factors. 

Simulations of such detailed protein models are extremely difficult. This is because the various competing inter- 
actions lead to multitude of local energy minima separated by high barriers. Hence, in the low-temperature region, 
canonical Monte Carlo or molecular dynamics simulations will get trapped in one of these minima and not thermalize 
within the available CPU time. Only with the introduction of new and sophisticated algorithms such as generalized- 
ensemble techniques [14], is it possible to alleviate this problem in protein simulations [15]. For this reason, our 
investigations rely on the use of one of these techniques, multicanonical sampling [16], where conformations with 
energy E are assigned a weight Wmu{E) oc l/n{E) {n{E) is that is the density of states). A simulation with this 
weight will generate a ID random walk in the energy space and lead to a uniform distribution of energy: 

Pmu{E) oc n{E) Wmu{E) = const . (7) 

Since a large range of energies are sampled, one can use the reweighting techniques [17] to calculate thermodynamic 
quantities over a wide range of temperatures T by 

- Jdxw-^E{x)) e-m-) ' 
where x stands for configurations and /3 = l/fcsT is the inverse temperature. Estimators for the multicanonical 
weights w{E) = n~^{E) = ex.p{—S{E)) can be calculated with the iterative procedures described in Refs. [6]. 



In our case we needed between 100,000 and 200,000 sweeps for the weight factor calculations. All thermodynamic 
quantities are then estimated from one production run of 2, 000, 000 Monte Carlo sweeps that followed 10, 000 sweeps 
for "thermalization" . Our simulations start from completely random initial conformations (Hot Start) and one Monte 
Carlo sweep updates every torsion angle of the peptide once. At the end of every 4th sweep, we store the total energy 
EToti the ECEPP/3 energy Eecepp/3, its partial terms Ec, Elj, Ehb and Etor, the solvation energy Esoiv, the 
corresponding end-to-end distance de-e, and the number nn {ns) of helical (sheet) residues. Here, we follow previous 
work [3] and consider a residue as helical if its backbone angle {(t>,tlj) are within the range (—70° ± 30°, —37° ± 30°). 
Similar, a residue is assumed to be "sheet-like" if {(t),il>) are within the range (-140° ±40°, 140° ±40°). 

III. RESULTS AND DISCUSSION 

We start with presenting our results for a single EKAYLRT molecule that is not interacting with other molecules. 
We display for this peptide in Fig. 1 the average helicity < nn > (T) as a function of temperature. Shown are 
data obtained in gas-phase (GP) and such for simulations that rely on a solvent accessible surface area term (ASA) 
to approximate protein- water interactions. We observe in both cases a steep helix-coil transition that separates a 
high-temperature region with little helicity from a low-temperature region where most of the residues are part of an 
a-helix. The location of this transition can be determined from the corresponding peaks in the specific heat C(T) that 
are drawn in the inlet. We find as the helix-coil transition temperature of EKAYLRT in gas phase T^f = 445 ± 15 K. 
The more pronounced peak for the solvated molecule indicates a temperature T^f^ = 340 ± 10 K that is considerably 
lower than the one in gas phase. Unphysiologically high helix-coil transition temperatures in gas phase, and their 
shift toward a more sensible temperature range when an implicit solvent is introduced, have been also observed in 
our earlier work on homopolymers [6,7]. 

We show in Fig. 2 as an example for the helical configurations that dominate below Thc the lowest energy config- 
uration found in a simulation of the solvated peptide {ETot = —69.6 kcal/mol). The lowest energy configuration in 
gas phase {Exot = Eecepp/3 — —28.0 kcal/mol) is also a helix (structure not shown). The energy of these helical 
structures is by » 25 kcal/mol lower than the ones of the lowest found "sheet-like" configurations: Exot = —43.8 
kcal/mol for the solvated peptide and Exot = Eecepp/3 = —3.1 kcal/mol for EKAYLRT in gas phase. 

The preference for helical structures can be also seen in Fig. 3a where we display the free energy AG at T = 300 
K as a function of helicity hh and "sheetness" hb- Note that for convenience we have chosen a normalization where 
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the minimum in free energy takes a value of zero. Both in gas phase and for the solvated molecule a funnel-like 
free-energy landscape is formed, with the free energy minimum at Uh = 5, i.e. for maximal helicity (since the two 
terminal ends are flexible and will usually not be part of an helix, a fully formed helix has a length uh = 5 instead 
of hh = 7). The absolute value of the free energy difference between coil and helix is much larger for the peptide in 
gas phase (AG ~ —5 kcal/mol) than it is for the solvated molecule (AG « —2 kcal/mol) indicating that the helix-coil 
transition is stronger for EKAYLRT in gas phase than for the molecule in an implicit solvent. This is in agreement 
with earlier work where we have found similar results for polyalanine chains [7] . The corresponding projection of the 
free energy landscape on the "sheetness" ns in Fig. 3b shows the opposite picture; the free energy increase with the 
number of residues whose backbone dihedral angles takes values that are common in a beta-sheet. Coil structures are 
at T = 300 K favored over sheet-like structures by AG « 5 kcal/mol in the implicit solvent and by AG « 8 kcal/mol 
in gas phase. 

The observed form of the free-energy landscape is caused solely by the intra-molecular interactions. This can be 
seen in Fig. 4 where we plot for solvated EKAYLRT the total energy Exot, the internal energy Eecepp/s and the 
solvation energy Esoiv as a function of temperature. Here, we have normalized all energy terms in such way that their 
value at nn — mO {ub = 0) is zero. Both ETot and Eecepp/z decrease with growing number of residues that are part 
of a helix while Esoiv increases (Fig. 4a). Hence, the protein- water interaction term opposes helix- formation. This 
result is reasonable as the protein-water hydrogen bonds compete with the intra-molecular hydrogen bonding in an 
a-helix and therefore weaken helix-formation in solution. However, the loss in solvation energy of A « 4 kcal/mol with 
helix-formation is small when compared with the gain in Eecepp/3 ~ —16 kcal/mol, and on average, a completely 
formed helix (hh = 5) has a total energy that is by AExot ~ —12 kcal/mol lower than a coil configuration (hh = 0). 
Not surprisingly, we observe the opposite behavior in Fig. 4b where we plot the same three energies as a function of 
"sheetness" Ub- Sheet-like configurations with large numbers Hb have higher internal energy E^cepp/z than such 
with ns = while the solvation energy Esoiv is lower. 

Hence, while at T = 300 K the protein-water interaction seems to favor strands and opposes helix-formation, the 
physics of our molecule is dominated by the intra-molecular energies that lead to a strong preference for a-helix 
formation. Fig. 5 indicates that this behavior is mainly due to the van der Waals interaction between the atoms 
in the peptide. In this figure, we display as a function of temperature besides the van der Waals term < Eydw > 
also the other partial energies that together make up Eecepp/z)'- the average electrostatic energy < Ec >, the 
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hydrogen-bond energy < Ehb > and the torsion energy < Etot >• 

Our results so far indicate that the peptide EKAYLRT has a intrinsic tendency to form hehces. Strands have of 
order w 30 kcal/mol higher free energies and are rarely observed. This result is independent on whether the molecule 
is in gas phase or simulated with an implicit solvent. However, EKAYLRT appears within proteins both in helices 
and /3-sheets. It follows that sheet formation has to be due to the interaction of the peptide with its surrounding. We 
conjecture that EKAYLRT forms a /3-sheet if it is in the proximity of another strand. Especially, we assume that this 
process also happens if the peptide is close to another EKAYLRT peptide that is already in a strand configuration. 
Unfortunately, the present version of SMMP does not allow the simulation of two interacting proteins. Hence, in 
order to test our conjecture, we have studied instead the peptide NH2-EKAYLRT-GGGG-EKAYLRT-COOH with 
the C-terminal EKAYLRT residues kept as /3-strand. The four glycine residues form a flexible chain that hold the 
two EKAYLRT-units together but allows their relative positions to vary. We refer to the so constructed peptide as 
molecule 'A'. 

The end-to-end distance dg-e is a measure for the separation of the two EKAYLRT chains. Our conjecture implies 
that for large values of de-e the N-terminal EKAYLRT assumes an a-helix while for small values of de-e (i-e. close 
proximity to the C-terminal EKAYLRT that forms a strand) it should assume a /3-sheet-configuration. We have 
therefore calculated from the multicanonical simulation of molecule 'A' the helicity and sheetness of the N- Terminal 
EKAYLRT at T = 300K. Both quantities are displayed in Fig. 6. Two regions are observed. For de-e >~ 16 A 
the N-tcrminal EKAYLRT chain forms a complete helix and strands are rarely observed. Hence, for these distances 
the N-terminal chain has a similar behavior as the isolated EKAYLRT-peptide. However, for decreasing end-to-end 
distance, the helicity also decreases and vanishes for de-e <~ 10 A. At the same time, the sheetness increases and the 
peptide forms a /3-shcct for dg-e ~ 5 — 6 A. Note that the average potential energy of helical configurations is with 
< Exot >= —24.9(1.6) kcal/mol within the errorbars equal to that of sheet-like configurations (< ETot >= —23.4(2.9) 
kcal/mol). 

In Fig. 7, the projection of the free-energy landscape at room temperature (T = 300 K) on the helicity and sheetness 
of the N-terminal EKAYLRT residues is drawn. For convenience, we have set in this figure the lowest-found value 
of the free-energy to zero as energies are only defined up to an additive constant. The contour lines are spaced by 2 
kcal/mol. The free-energy landscape is only plotted for values of G < 25 kcal/mol as values of the free energy grow 
rapidly outside of the drawn area. We observe again two minima, corresponding to fully formed helix and /3-strands. 
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Examples of configurations that correspond to the two minima are shown in Fig. 8. Both minima have comparable 
free energies and are separated by barriers of only 2 kcal/mol allowing an easy interchange between the two forms. 

In order to understand in more detail why EKAYLRT forms a /3-strand when close to a molecule that is already 
in a /3-sheet form, we have performed further simulations of NH2-EKAYLRT-GGGG-EKAYLRT-COOH holding now 
not only the C-terminal EKAYLRT-residues as a /3-strand but forcing also the four connecting glycine residues 
into a turn. We refer to the so defined peptide as molecule 'B'. The N-terminal EKAYLRT-residues are now by 
construction in close proximity to the C-terminal EKAYLRT-strand. Hence, we expect that at room temperature 
the N-terminal EKAYLRT chain will also form a /3-strand. This conjecture is supported by Fig. 9 where we plot the 
average "sheetness" ris of the N-terminal EKAYLRT-residues as a function of temperature. Both in gas-phase and 
for simulations with a solvent accessible surface term, we find that on average more than 5 of the 7 residues are part 
of a sheet-like structure. < ns > decreases smoothly with growing temperature and the maximum in the specific heat 
is shallow. The transition is more pronounced for the peptide in an implicit solvent than for the one in gas-phase, 
and shifted toward lower temperatures. 

Unlike our previous simulations where the glycine residues could freely move, a large percentage of configurations 
are now at room temperature in a sheet-form. The increased statistics of these configurations allows for a better 
analysis of the factors that help to overcome the intrinsic propensity of EKAYLRT to form an a-helix and lead to a 
/3-sheet. Table 1 lists the diff'erences of various energies between structures where the N-terminal EKAYLRT unit is a 
/3-strand with structures where these residues form an a-hclix. Values are listed for the whole molecule 'B' and such 
restricted to the N-terminal EKAYLRT chains. Also listed are the differences of both terms. The latter quantity is a 
measure for the interactions between these seven residues and the rest of the molecule (that is kept fixed). 

We see from this tabic that at T = 300 K configurations with the N-tcrmiiial EKAYLRT chain in a sheet are 
energetically favored by 7 kcal/mol over such where these residues form an a-helix. This energy-bias is found for all 
partial energies with the exception of the solvation energy term Esoiv and the torsion energy term Etot- While their 
values seems to indicate a slight preference for helical over sheet-like configurations, they are within the errorbars 
compatible with zero suggesting that both terms show no preference for one of the two forms. 

Because of the bias in the internal energies _E'^o'«cu;e 's q£ ^j^g -^i^ole molecule 'B', /3-sheet-like conformations of 
the N-terminal EKAYLRT chain dominate at room temperature. However, the behavior of the various energy terms is 
different when one considers only the contributions by these seven residues. With the exception of the solvation energy 
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Eg^^'^^^^ , that favors /3-strands, all energy terms favor now a helix. On average, helical structures have at T = 300 
K a 14.5(1.4) kcal/mol lower energy e^kaylrt ^^jjg^n strands when only the interaction between atoms in these seven 
residues are considered. Hence, their behavior is qualitatively the same as for the isolated EKAYLRT-peptide where 
we also observed a strong bias toward helical conformations. Again, we find also that the van der Waals energy E^dw 
is the dominant term. It follows that the /3-sheet configurations that dominate when the EKAYLRT residues are build 
into molecule 'B' are caused by the interaction between this chain and the "background" of the rest of the molecule. 
Since energies are additive, we can calculate this "background" -field by 

j^Background j^Molecule^ B' jjjEKAYLRT ^gj 

The strength of the interaction between the peptide and the "background field" of the rest of molecule 'B' can be seen 
from the large energy- difference of /^^j^^^^'^sround _ _21.8(1.6) kcal/mol by that these interactions favor a strand. 
The main contribution comes from the van der Waals term (^AE^^^^^^ound ^ 15.7(1.1) kcal/mol) which is almost 
three times as large as the electrostatic and torsion energy terms. Note that Eq. 9 tell us also that the "background" 
given by the fixed parts of molecule'B' raises the solvation energy difference AEg^^^^^'^ = —3.1(2) kcal/mol of the 
EKAYLRT-chain by AE^f^''®™""'^ = 3-7(3) kcal/mol to a value of AE^f^^"''^^''^' = 0.6(3) kcal/mol for the whole 
system. This is because the term AEg^^^^^'^ is due to the competition between hydrogen-bond formation in an 
a-helix and hydrogen-bond formation between the peptide and the surrounding water. However, in Molecule 'B' the 
peptide is geometrically constraint in such a way that this competition is replaced by one between hydrogen-bond 
formation in an a-helix of the N-terminal EKAYLRT on one side, and formations of hydrogen-bonds between the 
peptide and the C-terminal EKAYLRT-residues on the other side (see also the opposite sign in the terms AE^^^^^^^ 
and A£;f ^''^''"""'^ in table 1). As a result, both the solvation energy difference AEf^l^"'^^'' and hydrogen-bond 
energy difference AE^g""''^'' are marginal. Instead, the preference of /3-sheet configurations for EKAYLRT in 
the "background" of the fixed rest of Molecule 'B' seems to be mainly due to the electrostatic and van der Waal's 
energies. This is reasonable: a /3-sheet conformation allows for an average closer distance between the atoms of the 
N-terminal EKAYLRT-chain and the existing /3-strand of the C-terminal EKAYLRT-residues, decreasing in this way 
the van der Waals energy. At the same time, the alignment of the two /3-strands leads also to a favorable alignment 
of the dipole moments associated with each residue lowering therefore the electrostatic energy. We conjecture that 
without the stereometric constraints imposed by the connecting glycine residues the two strands would move together 
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and aggregate as the energy gain increases with decreasing distance between them. 

Our above presented results for Molecule 'A' and Molecule 'B' suggest auto-catalytic properties for EKAYLRT: if 
the peptide forms a strand, in becomes energetically favorable for other nearby EKAYLRT molecules to transform 
themselves into a sheet (instead of the normally preferred helix), and eventually to aggregate with the first one. This 
behavior is similar to the mechanism thought to be responsible for the outbreak of neuro-degenerative illnesses such 
as Alzheimer's or the Prion diseases. Outbreak of theses illnesses is associated with the appearance of a miss-folded 
structure that differs from the correctly folded one by a /J-sheet instead of an a-helix. The miss-folded structure is 
thought to be auto-catalytic, that is its presence leads to a structural transition by which the correctly folded (helical) 
structure changes into the harmful /3-sheet form. Hence, peptides that contain the sequence of amino acids EKAYLRT 
can serve as simple models to study these a — > /3-transition and the mechanism of Prion diseases. For instance, our 
investigation suggest that the formation of /3-sheets can be minimized by shielding the surface area of already existing 
/3-sheet forms minimizing in this way the van der Waals interaction. Another possibility may be to introduce metal 
ions that alter the electrostatic interaction decreasing in this way the energy bias toward /3-sheets. 

IV. CONCLUSION 

We have performed multicanonical simulations of peptides that contain the sequence of amino acids EKAYLRT. We 
find that the EKAYLRT-peptide itself has both in gas-phase and in solution an intrinsic tendency to form an a-helix. 
However, the peptide assumes a /3-sheet form when close to another strand. The transition from an a-helix toward a 
/3-sheet is caused by strong van der Waals und electrostatic energy terms that favor the /3-sheet form when EKAYLRT 
is in close proximity to another strand. This "auto-catalytic" property of EKAYLRT, that induces strand-formation 
in other EKAYLRT molecules when in a /3-sheet configuration, suggests that the EKAYLRT based peptides can serve 
as a simple model for the a — > /3-transitions and successive aggregation that are supposed to be related to the outbreak 
of various illnesses such as Alzheimer's or the Prion diseases. 
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Tables: 





Aloloculo -B' 


X-tcrmiual EKAVLllT residues only 


Background 


AExot 


-7.2(9) 


14.5(1.4) 


-21.8(1.6) 


AEsoiv 


0.6(3) 


-3.1(2) 


3.7(3) 


AEel 


-4.3(3) 


0.6(1) 


-4.9(3) 


AEydW 


-3.2(8) 


12.5(9) 


-15.7(1.1) 


AEhb 


-0.8(2) 


4.0(2) 


-4.8(3) 


AEtot 


0.5(4) 


0.5(4) 


0.0 



TABLE I. Energy differences between "sheet" and helix configurations at room temperature for various energy terms as 
calculated from a multicanonical simulation of molecule 'B'. 
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Figure Captions: 



Fig. 1 The average number < uh > of helical residues as a function of temperature T for EKAYLRT in gas phase 
(GP) and simulated with an implicit solvent term (ASA). The specific heat C{T) as function of temperature T 
is displayed in the inlet. All results rely on multicanonical simulations of 2,000,000 sweeps each. 

Fig. 2 Lowest energy configuration of EKAYLRT as found in a multicanonical simulation of 2,000,000 sweeps using 
an implicit solvent to approximate the peptide-water interactions. 

Fig. 3 The free energy AG at T = 300 K as a function of {bottom) helicity and {top) sheetness ub for EKAYLRT 
in gas phase (□) and simulated with an implicit solvent term (A). The free energy is normalized in such a way 
that its minimum value is set to zero. All results are calculated from a multicanonical simulation of 2,000,000 
sweeps. 

Fig. 4 The average total energy < -Eto* > (1^)) intra-molecular energy < Eecepp/3 > (^) and solvation energy 
< Esoiv > (°) of EKAYLRT at T = 300 K as a function of {bottom) helicity uh and {top) sheetness ns- AH 
energies are normalized in such way that their value at hh = (ns = 0) is zero. All results are calculated from 
a multicanonical simulation of 2,000,000 sweeps using an implicit solvent model to approximate peptide-water 
interactions. 

Fig. 5 The average partial energies < Ec >, E^iiw^Ehb and Exor that together make up the ECEPP/3 energy 
Ebcepp/3 as a function of temperature T. All terms are normalized in such way that their value for T = 1000 
K is zero. All results are calculated from a multicanonical simulation of 2,000,000 sweeps using an implicit 
solvent model to approximate peptide-water interactions. 

Fig. 6 The average helicity < nn > and sheetness < ub > &t T = 300 K of the N-Terminal EKAYLRT residues 
as a function of the end-to-end distance dg-e- AH results are calculated from a multicanonical simulation of 
2,000,000 sweeps. 

Fig. 7 The free-energy landscape of molecule 'A' at room temperature {T = 300 K) as a function of helicity < uh > 
and sheetness Kub >■ The global minimum is set to zero and the contour lines are spaced by 2 kcal/mol. 
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Fig. 8 Low-energy configurations of molecule 'A' as found in a multicanonical simulation of 2,000,000 sweeps. The 
one in (a) is the lowest-energy configuration where the N-terminal EKAYLRT-residues form an a-helix; the one 

in (b) where they form a /3-sheet. 

Fig. 9 The average "sheetness" < > of the N-terminal EKAYLRT-residues of molecule 'B' as a function of 
temperature. 
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